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Abstract 

We suggest the method for group classification of evolution equations admitting 
, nonlocal symmetries which are associated with a given evolution equation possessing 

nontrivial Lie symmetry. We apply this method to second-order evolution equations 
in one spatial variable invariant under Lie algebras of the dimension up to three. 
Ch . As a result, we construct the broad families of new nonlinear evolution equations 

possessing nonlocal symmetries which in principle cannot be obtained within the 
classical Lie approach. 
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i> ■ 1 Introduction 

t— i ■ 

Classical Lie symmetries has become an inseparable part of the modern mathematical 
' physics toolkit for analysis of linear and nonlinear differential equations pQ-jO]- However, 

with all its power, versatility and universality the Lie symmetry approach has its limits, 
which prompted numerous researchers to look for possible generalizations of the concept 
of classical symmetry in order to be able to handle differential equations that do not 
possess Lie symmetries. 

One of the promising generalizations of the concept of the classical symmetry is al- 
lowing for a generalized symmetry to depend on integrals of dependent variables, in con- 
trast to Lie symmetry which can only involve independent, dependent variables and their 
derivatives. In this way nonlocal symmetries have been introduced into the modern math- 
ematical physics. 

There was significant progress in studying nonlocal symmetries of linear differential 
equations (see, e.g., [7] and the references therein). Much less is known about nonlocal 
symmetries of nonlinear differential equations. 

A possible way to approach the problem of constructing nonlocal symmetries of partial 
differential equations is utilization of the concept of quasi-local symmetry introduced in 
[8], [9] . Quasi-local symmetry is the one obtained from a classical Lie symmetry through 
nonlocal transformation. Another approach is based on the concept of potential symmetry 
[10]- [12] . However, as we established recently, potential symmetries of evolution type 
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equations are quasi-local ones in a sense that there is a nonlocal transformation reducing 
those to point or contact symmetries [T3] . 

The principal motivation for this paper was a need for a more systematic approach 
to classification of partial differential equations admitting quasi-local symmetries. It is 
natural to expect that such an approach should rely on Lie symmetries of equations under 
study We demonstrate in the paper that it is indeed possible to develop the regular 
and systematic classification scheme using the results of group classification of nonlinear 
second-order equations [HJ [15] and the ideas suggested in [T3l [16] . 

We utilize the results of group classification of the general second-order evolution 
equation in one spatial variable 

u t = f{t,x,u,u x ,u xx ) (1) 

to construct nonlocal symmetries of the associated evolution equations. In a what follows, 
we denote the class of partial differential equations (p]) as E. 

The paper is organized as follows. In Section 2 we give the necessary definitions and 
present the detailed description of our approach to classifying nonlocal symmetries of ([I]). 
Section 3 is devoted to application of the approach to equations from E invariant under 
the Lie algebras of the dimension s < 3. The last section contains brief discussion of the 
obtained results. 



2 Description of the method 

We begin by introducing the differential transformation which plays a key role in our 
subsequent considerations 

CP(t, x, u) = (i, x, u) = (t, x, u x ). (2) 

Note that the inverse of differential transformation ([2]) reads as 

t — t , x = x, u = u. (3) 

Hereafter d~ x is the inverse of the differentiation operator d x , i.e., d x <9~ x = d^ 1 d x = 1. 
Consider a subclass of evolution equations of the form ([1]) 

u t = f(t,x,u x ,u xx ) (4) 

admitting the group of translations by u. 

Differentiating with respect to x and applying ^3 yields 

ui = d x f(t,x,u,u x ). (5) 

Consequently, ^3 maps a subclass of equations (jlj) from E into another subclass of equa- 
tions (jS]) from E. Remarkably, this very transformation is responsible for all potential 
symmetries of equations (pP) [13] . 
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Indeed, suppose that PJ admits Lie transformation group 



t' = T(t,6), x' = X(t,x,u,6), u' = U(t,x,u,6), (6) 

where 9 is the group parameter. Note that is the most general one-parameter Lie group 
that can be admitted by an equation from the class E (see, e.g., Computing the 

first prolongation of formulas ([6]) we derive the transformation rule for the first derivative 
of u 

du' _ U u u x + U x 
dx' X u u x + X x ' 

Now the image of transformation group ((5]) under the action of takes the form 

? = T(i, 8), x' = X(t,x,u,9)), u'= ^ u _ + p (7) 

X u u + X x 

with u = d^u and U = U(i,x,u,9). Consequently, if the right-hand sides of ([7]) contain 
nonlocal variable u, then ([7]) is nonlocal symmetry of (JH). According to [13] . transforma- 
tions ([7]) include the variable u if and only if X u ^ or U? u + U 2 U ^ 0. 

Summing up we conclude that provided symmetry group ([6]) of equation (jlj) obeys one 
of the constraints 

X u ± or U 2 XU + U 2 UU ± 0, (8) 

then *P maps this symmetry into nonlocal symmetry (I7J) of equation ([5]). This nonlocal 
symmetry is called the potential symmetry of (J5]). 

Constraints ([8]) can be written in terms of the coefficients of the infinitesimal operator 
of group §6§ 

S = r(t)d t + f (i, x, u)d x + rj(t, x, u)d u . (9) 

If 

£^0 or t^ + ^^O, (10) 

then symmetry S 1 is mapped by ^3 into nonlocal symmetry |13j . 

Despite the fact that equation (jl]) is rather particular case of the general evolution 
equation ([T|), the scheme outlined above provides a very general framework for group 
classification of nonlocal symmetries. 

Let the symbol S stand for the infinite-dimensional Lie algebra of differential opera- 
tors ([H]) and the symbol I stand for the infinite-dimensional Lie algebra spanned by the 
operators of the form 

Q = £(t,x,u)d x + r](t,x,u)d u . (11) 

As direct computations shows, [X, X] C X and [X, S] C X, which means that X is the ideal 
in the Lie algebra S. 

If an equation of the form ([T]) possesses Lie symmetry from X, then it can always be 
reduced to the form (j4]) by a suitable local transformation of variables [2]. Consequently, 
the procedure for classification of nonlocal symmetries described above is applicable to 
any equation from the class E admitting at least one symmetry from X. 
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In the paper [T3] we suggest a generalization of the potential symmetry approach 
which allows for the use of more general differential transformations to construct nonlocal 
symmetries. Our generalization is based on the idea put forward by Sokolov [16J. He 
suggests that a differential transformation is to be looked for in the form 

D(t,x,u) = (i,x,u) = (t,ui(t,x,u,u x ,u xx , . . .),u 2 (t,x,u,u x ,u xx , ...)), (12) 

Ui and uj 2 being differential invariants of the symmetry group of the evolution equation 
under study. For the case when D = we have u>i = x and u 2 = u x , x and u x 
being the simplest differential invariants of the group of displacements by u generated 
by the infinitesimal operator d u . Evidently, evolution equation (J4]) admits symmetry d u . 
Sokolov proves in [TH] that for any contact symmetry of evolution equation of the form 
([T]) preserving the temporal variable t the corresponding differential transformation maps 
equation in question into E. 

Now we are ready to formulate the group approach to classifying nonlocal symmetries 
of equations from the class E. 

Let equation (TjTJ be invariant under the r-dimensional Lie algebra C r = (ei, . . . ,e r ). 
Next, suppose that C r contains a subalgebra l s = (ei,...,e s ), s < r such that T s C 
X. Saying it another way, the basis elements of the algebra X s have the form (fill) . Let 
t,u>i,uj2 be functionally independent differential invariants of X s . 

Making differential transformation (fl2l) with so chosen U)i,u)2 we map the invariant 
equation ([T]) into another equation from E. If the transformation group generated by 
the operators e s+ i,...,e r obeys one of the inequalities (1121) . then its image under the 
differential transformation D is a nonlocal symmetry of the transformed equation. 

Direct application of the transformation of the form (JT]), especially for the case when 
Ui,U2 depend on higher derivatives, may become quite challenging computational prob- 
lem by itself. That is why we suggest a modification of the approach above, which is to 
represent D as a superposition of elementary differential transformations. These elemen- 
tary transformations are obtained as a superposition of a local transformation reducing 
the individual basis operators of the Lie algebra X s to the form d u and of the differential 
transformation *}1 

We begin with the local change of variables, %\, reducing the first basis operator e% 
to the form d u and after that we apply the transformation ^3. Next, we recompute the 
remaining basis operators in the new variables t, x, u x . 

To construct the image of infinitesimal operator (ITTj) under the transformation *p we 
compute the coefficient of d Ux in its first prolongation 

C(t, X, U, U x ) =T] X + U X (T] U - £ x ) - u 2 x ^ uu . 

Next, we replace u with d~ l u and u x with u so that the transformed infinitesimal operator, 
Q, takes the form 

Q = x, u)d x + ((t, x, d^u, u)d u . (13) 

If there is a basis element, e^, which does not depend on the nonlocal variable d^u, then 
we can repeat the procedure and reduce e« to the form d u and then apply ^3 and so on. 
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The above described procedure has two possible outcomes. The first one is when we are 
able to perform all s steps and get an evolution equation which is the image of the initial 
X s -invariant equation under transformation ffT2l) . The resulting differential transformation 
takes the form 

1 

%=H^o%) = <poT s o<poT s _ 1 o...o<poX 1 . (14) 

i=s 

Applying X to JT]) yields the sequence of s evolution equations from the class E together 
with the sequence of their symmetries, some of which might be nonlocal. 

The second outcome is when at some step p < s all the remaining basis operators 
contain the nonlocal variable d~ l u and we cannot proceed any further. In this case we 
obtain the sequence of transformations 

1 

1=11(^0%) = ^o1 p o^o% p _ 1 o...o^o% 1 (15) 

i=p 

that maps the initial X s -invariant evolution equation into equation, which admits at least 
s — p nonlocal symmetries of the form (fT5|) . 

We summarize the above speculations in the form of the algorithm of classification of 
evolution equations possessing nonlocal symmetries. 

The starting point of our classification algorithm is evolution equation ([1]) invariant 
under r-dimensional Lie algebra C r C S. Suppose also that the algebra C r contains a 
subalgebra 1 S C X of the dimension 1 < s < r. Now to classify evolution equations 
admitting nonlocal symmetries one needs to 

• classify X s -inequivalent subalgebras J' 1 , . . . , J N of the algebra X s , 

• compute the sequence of transformations f[T41) mapping the initial invariant equation 
into E for each subalgebra J 1 , i = 1, . . . , N and generate the corresponding sequence 
of evolution equations and their invariance algebras, 

• verify for each invariance algebra whether its basis elements satisfy one of the in- 
equalities ffTU]) . If this the case, then the remaining elements of the sequence of 
evolution equations possess nonlocal symmetry. 

Upon completing the three steps of the algorithm above ones obtains the set of evo- 
lution equations ([T|) that admit nonlocal symmetries. 

3 Applications 

In [T4"l [To*l [T7] we obtain exhaustive description of inequivalent Lie algebras, which are 
symmetry algebras of equations of the form ([1]). This fact facilitates application of the 
algorithm presented in the previous section in its full generality. However, in this paper we 
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restrict our considerations to subalgebras of the symmetry algebra of (pQ) of the dimension 
s < 3. 

The list of inequivalent realizations of one-dimensional symmetry algebras of equation 
(JTJ) consists of one representative 

Ai = (d u ) : u t = fit, x, u x , u xx ). 

This case has already been considered in the previous section, the algebra A\ giving rise 
to the transformation 

According to [17] the class of operators, X, contains four inequivalent realizations of 
two-dimensional symmetry algebras of ([T]). These realizations are listed below together 
with the corresponding invariant equations 

A\ = (d x , d u ) : ut = f{t, u x , u xx ), 

= {-^9 U ,0 U ) '. Ut — f(t,X,U xx ), 

A2 = ( xO x uO u , U ) . Ut — x f (t, u x , xu xx ) , 
A\ = (-ud u , d u ) : u t = u x f(t, x, u~ l u xx ). 

We start by analyzing differential transformations associated with the algebra A\. 
Applying the transformation *P to the equation u t = f(t, u x , u xx ) we get 

where / = f(t, u, u x ). Note that we dropped the bars. 

The second basis operator of A\ reads now as d x . To reduce it to the the form d u we 
apply the hodograph transformation 

%i(t,x,u) = (i,x,u) = (t,u,x). (16) 

The transformed equations has the form 



f' , U XX r 

u t = ~J Z1 H Jz 2 , 



where / = f(t,zi,Z2) = f(t,u, —u x l ). Applying the transformation yields the final 
form of the transformed evolution equation 



d 



(-U + -f z 

V u 



OX V U 

Here / = f(t, z h z 2 ) = f(t, u, -u' 1 ). 

Turn now to the algebra A\. Applying the transformation to the invariant equation 
u t = f(t,x,u xx ) yields 

Ut fx U xx f Ux 

with / = f(t,x,u x ). As usual we drop the bars. The second basis element of A\ takes 
the form d u and we can directly apply the transformation ^3 thus getting 

Ut fxx 2>U x f xu ~\- U x f uu ~\- fuU xx - 
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* dx V 



Consider next the algebra A\. Applying the transformation *p to the equation u t = 

xf(t,u x ,xu xx ) we have 

u t = f + xu x f zi + x(u x + xu xx )f Z2 (17) 

with / = f(t,Zi,z 2 ) = f(t,u,xu x ). The second basis operator of the algebra A\ reads 
now as xd x . Applying the hodograph type local transformation 

1 2 (t,x,u) = (i,x,u) = (t, u, Inx) (18) 

we reduce it to the form d u . Transforming equation (fTTj) accordingly we get 

Uf "icfzi (pC U x ^ra %U X •^)y^2> 

where / = f(t, zi, z 2 ) = f(t, x, u' 1 ). Now we can apply the transformation thus getting 

— f-xfzj. + (x 2 u~ 2 u x + xu' 1 - x)f Z2 ^j . 

with / = f(t,z 1 ,z 2 ) = fit.x.u^ 1 ). 

Finally, we turn to the algebra A\. Applying the transformation to the invariant 
equation u t = u x f(t, x, u^u^) yields 

u t = u x f + uf zi + {u xx - u~ x u 2 x )f Z2 , (19) 

where / = f(t, Z\, z 2 ) = f(t, x, u~ lr u x ). The second basis element of the algebra A\ is still 
of the form ud u , so we need to reduce it to the form d u . To this end, we perform the local 
transformation 

%z{t, x,u) = (t,x,u) = (t,x,\nu). (20) 
As a result, (fT9~j) takes the form 

Ut — u x f + f Zl + u xx f Z2 
with / = f(t,zi,z 2 ) = f(t,x,u x ). Applying the transformation *P finally yields 

Ut u x f -\- uf x -\- uu x f u -\- f xx -\- 2u x f xu -\- u x f uu -\- u xx f u , 
where / = f(t, x, u). 

Summing up we present the four sequences of differential transformations that map 
four subclasses of the class of second-order evolution equations E into E. 

• Differential transformation *P o % 1 o <p maps equation 

u t = f(t,u x ,u xx ) 

into 

d ( u x 
Ut -dx~\ U + V /z2 
with / = f(t,z u z 2 ) = f(t,u, -l/u). 
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Differential transformation *p o ^3 maps equation 

Ut = f(t, X, u xx ) 

into 

Uf fxx 2'U x f xu U x f U u -\- fuUa 

with / = f(t,x,u). 

Differential transformation *p o % 2 ° maps equation 

u t = xf(t , U x , XU XX ) 

into 



d ( 

Ut= dx \ x f Zl + ( x2u %Ux + xu ~ x ^f- 



with / = f(t, z x ,z 2 ) = /(t, x, u 1 ). 
• Differential transformation o % 3 o ^3 maps equation 

«t = U x f{t,X,U~ X U xx ) 

into 

^ — w x / + It/j; + uu x f u + / zx + 2u x f xu + u 2 x f uu + u xx f u , f = fit, x, u). 

In the above formulas the transformations Xi,^, X3 are given by f|T6|) . (IT7j) and (|20p . 
respectively. 

The class of first-order differential operators I contains twenty two inequivalent real- 
izations of three-dimensional symmetry algebras of ([I]) [17]. These realizations are listed 
below together with the corresponding invariant equations. 
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— xd x — ud u , d u , xtd x ) : u t = —xt 1 u x In u x + xu x f(t, xu x u a 
—xd x — ud u , d u , xd x ) : Ut = xu x f(t, xu x l u xx ), 
-xd x - ud u , d u , xd u ) : u t = xf(t, xu xx ), 
-ud u , d u , d x ) : u t = u x f(t, u x l u xx ), 

d u , d x , xd u + td x ) : u t = -^ u l + /(*> u xx ), 

d u , d x , xd u ) '■ Ut = f(t,u xx ), 
xd u , d u , x 2 d x + xud u ) ■ u t = xf(t, x 3 u xx ), 

d x , d u , (x + u)d x + ud u ) ■ u t = u x exp(u~ 1 )/ (t, u x 3 u xx exp(u~ v 

xd u ,d u ,x 2 d x + u(l + x)d u ) : u t = xexpf^-aT 1 )/^, ot^ exp^ -1 ) J . 
(d x , d u , xd x + ud u ) ■ u t = u xx f(t, u x ), 
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■Af = ( 9 x, d u , xd x - ud u ) : u t = ul /2 f(t, u x 3/2 u xx ), 

Af = (xd u , d u , -2xd x - ud u ) : u t = x 1/2 f(t, x 3/2 u xx ), 

A 1 , 3 = (d x ,d u ,xd x + qud u ) : ut = ufi-Vf(t,u x 2 ^-Vu xx ), 0<\q\< 1, 

A\ 4 = (xd u ,d u ,x(q-l)d x + qud u ) : u t = x«^f(t, x^^u^), g^O,±l, 

Al 5 =(d x ,d u ,ud x -xd u ) : u t ={l + u 2 x ) 1 / 2 f(t,(l + ul)- 3 / 2 u x ^, 

Af = (xd u , d U7 (1 + x 2 )8 x + xud u ) : Ut=(l + x 2 ) 1/2 f (t, (1 + x 2 ) 3/2 u xx ^j , 
^4g 7 = (d x , d u , (qx + u)d x + (—x + qu)d u ) : u t — (1 + m^) 1//2 exp(— q arctanu x ) 

xf(t, (1 + u 2 x )~ 3/2 exp(-g arctanwa,)^^ , g > 0, 
Alj 8 = (xd u , d u , (1 + x 2 )^ + (x + q)ud u ) : u t = (1 + x 2 ) 1//2 exp(arctanx) 

xf(t, (1 + x 2 ) 3 / 2 exp(— g arctanx)^^ J , g ^ 0, 

^4g 9 = (<9 U , cos ud x + tan x sin ud u , — sin -uc^ + tan x cos ii<9„) : 

Ut = (sec 2 x + w 2 ) 1 ^ 2 /(t, (« M cosx — (2 + w 2 cos 2 x)^ sinx)(l + u 2 x cos 2 x)~ 3 / 2 ), 
.A 20 = (2ud u - xd x , -u 2 d u + xud x , d u ) : u% = xu x f(t, x~~ 5 u x 3 u xx + 2x~ 6 u~ 2 ), 
Af = (2ud u - xd x , (x" 4 - u 2 )d u + xud x , d u ) : u t = x~ 2 (4 + x 6 u 2 ) 1/2 

xf(t, (4 + x 6 u 2 x )- 3/2 (x 4 u xx + 5x 3 u x + ^x 9 m 3 )), 

.4 22 = (2ud u - xd x , -(x~ 4 + w 2 )<9 u + xud x , d u ) : ^ = x~ 2 (x 6 w 2 - 4) 1/2 

x/(t, (x 6 u 2 - 4)- 3 / 2 (xV c + Sx 3 ^ - ^x\ 3 )). 

Analysis of the sequences of differential transformations associated with the above algebras 
is similar to the case of two-dimensional algebras. We skip the derivation details and give 
the final result, the differential transformations (114)) . (1151) . 

For the algebras A\- A\, .A3 - A\ A The sequence of transformations X reads as 

T=^oT 2 ofpoT 1 o<p. (21) 

where 





%i(t, x, u) 


= {t,u, 


-In 


x), X 2 (t,x 


u) = (t, I«, 


t _1 lnx) 


i 


%i(t, x, u) 


= (t,u, 


-In 


x), X 2 (t,x 


■u) = (t, xu, 


lnx). 


A 3 


%i(t, x, u) 


(t, X, 


u), 


X 2 (t,x, w) = 


= (£, xm, lnx^ 




i 


%i(t, x, u) 


(t, X, 


-In 


u), X 2 (t,x 


w) = (t, xu, 


lnx). 


A 3 


%i(t, u, x) 


= 


x), 


X 2 (t,X,?i) = 


- (t, u, x). 






%i(t, x, u) 


= (t,x, 


u), 


X 2 (t,x, w) = 


= (t,u,x). 




AlO 


: % 1 (t,x,u] 


= (t,u 




X 2 (t,x, w) 


= (t, x, In u) 




A 3 




= (t,u 




X 2 (t,x,w) 


= (t, x 3 u 2 , - 


-lnw). 
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A?: 




u) 


= (t,u 




% 2 (t,x 


,u) 


= (t,x 3 u 2 , - 


1 

3 


\nu). 




A 13 ■ 


%x(t,x 


u) 


= (t,u 


,x), 


1 2 (t,x 


,u) 


= (t, x q -\ q - 


-l 


g-2 lnM) ' 






%l(t,X, 


u) 


= (t,u, 


x), 


% 2 (t,x, 


u) 


= (t, u, \nu), 




q = 2, 




Af: 


%i(t,x 


u) 


= (t,x 


,«), 


% 2 (t,x 


,u) 


= (t, x q -\ q - 


-l 


g-2 lnM) ' 






%l{t,X, 


u) 


= (t, X, 


u), 


% 2 {t,x, 


u) 


= (t, u, \nu), 




9 = 2, 





Now if an equation invariant under one of the algebras from the list ^.3 - A%, A\° - 
A\ A admits additional Lie symmetry satisfying ( flOj) . then every member of the sequence 
of equations obtained by successive application of transformations ^oljo^oljo^ 
possesses nonlocal symmetry Moreover, if we are able to construct an exact solution of 
some equation from the sequence, we can map this solution into solutions of the remaining 
equations from the sequence. 

For the algebras A\ - A\, Af - Af, and .4J - A% 2 the algorithm stops after the second 
iteration, since the third symmetry operator turns into a nonlocal one. Below we give the 
corresponding sequences of transformations. 



Al: 






At: 










Af 


: ^o% 1 


0% 


Af 


: %Soy, 






: (Pol! 


0% 


A? 


: ^0^, 




Af 


: ?oI 2 


0% 


A 21 


: ^oT 2 


0% 


Af 


: ?oI 2 





where ^ 

Ti(t, x, u) = (t, m, x), ^ 2 (t, x, u) = (t, x 3 u, - ln-u). 

o 

Applying the above sequences of transformations to the corresponding invariant equations 
yield evolution equations of the form ([1]) that admit nonlocal symmetries. The explicit 
form of those symmetries is very cumbersome, therefore we present the nonlocal sym- 
metries obtained through the action of ^5 only. In addition, we give the corresponding 
transformed equations. 

Al : e 3 = x 2 d x + (v- xu)d U} u t = ^- (xf(t, x 3 u x )) , 
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A 3 '. 


e 3 : 


v4 3 : 


e 3 ; 


A 15 ■ 


e 3 


A 16 ■ 


e 3 


A. 3 : 


e 3 



d_ 

dx 



= x 2 d x + ((l-x)u + v)d u , u t = -^(xexp(-x ^/(i, xexp(x V^)). 



9a; 



<9x V 
9 



x/(i, (1 + w 2 )~ 3/2 exp(-garctanw)w x )), (22) 

^ 8 : (1 + x 2 )d x + ((q - x)u + v)d u , u t = — ({1 + x 2 ) l/2 exp(arctanx) 

xf(t, (1 + x 2 ) 3 / 2 exp(—q aictaax)u x )^j , 

Al° : e 2 = xvd x — (xu + 3uv)d u , u t = — (xuf(t, x~ Fj u~ 3 u x + 2x~ 6 u~ 2 ) S J , 

d ( 

Af : e 2 = xud x - (4a;~ 5 + xu + 3uv)d u , u t = — ( x~ 2 (4 + iV) 1/2 
xf(t, (4 + x 6 u 2 )-^ 2 {x A u x + 5x 3 u + liV))j, 

^22 . ^ = + ( 4a ,-5 _ xu _ 3 uv )d u , U t = — I X~ 2 (X 6 U 2 - 4) 1/2 

x/ (i, (xV - 4)~ 3 / 2 (a;V E + 5x 3 u - ^V)) ^ . 

Here i> = d" 1 ^ is the nonlocal variable. 

Finally, when computing transformations associated with the algebra A!, 9 we have to 
stop after the first iteration, since the second and third basis operators are mapped by *}3 
into nonlocal symmetries 

e 2 = cos ud x + (u sin v + sec 2 x sin v + u cos v tan x)d u , 

e 3 = — sinud x + (u cosv + cos v sec 2 x — u sint> tanx)<9 u , (23) 

where v = d~ l u. So that the sequence of differential transformations boils down to a 
single transformation T = ^3. And what is more, the evolution equation invariant under 
the nonlocal symmetries e 2 and e 3 reads as 

df 

u t = (sec 2 a; tan a; + uu x ) (sec 2 x + u 2 )~ l/2 f + (sec 2 a; + u 2 ) 1/2 -^- 

with 

f — f(t, (u x cosx — (2 + -u 2 cos 2 x)-usina;)(l + u 2 cos 2 x)~ 3 ^ 2 j . 
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Let us emphasize that symmetries (1221) . (I2"3"j) contain nonlocal variable and, conse- 
quently, cannot be obtained within the framework of the traditional infinitesimal Lie 
method. 

4 Concluding Remarks 

We develop the regular method for group classification of evolution equations admitting 
nonlocal symmetries. The starting point of the method is an evolution equation that 
possesses nontrivial Lie symmetry. The source of nonlocal symmetries are classical sym- 
metries subjected to nonlocal transformations. 

One of the by-products of our approach are the sequences of evolution equations 
related to the initial invariant equation. So that if we succeed in constructing a (general 
or particular) solution of the sequence member then it is possible to map this solution 
into solutions of the remaining equations from the sequence. 

Of special interest is the case when one of the members of the sequence is the linear 
evolution equation, since this means that the remaining equations are linearizable. More- 
over, a direct application of the method to linear evolution equations admitting nontrivial 
Lie symmetries is also of importance, since it yields sequences of linearizable evolution 
equations. 

These and related problems are under study now and will be reported elsewhere. 
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